DIRICHLET PROBLEMS OF A QUASI-LINEAR ELLIPTIC 

SYSTEM 



GONGBO LI AND LI MA 

Abstract. We discuss the Dirichlet problem of the quasi-linear elliptic system 
-e- /(E/) &(e /((,) vf/)+ \J'{U)\X7U\ 2 =0, in O, 

u\ m = <p. 

Here Q. a smooth bounded domain in R n , f : R N — > R is a smooth function, 
U : CI — > R N is the unknown vector-valued function, 0:0—* i? w is a given vector- 
valued C 2 function, /' is the gradient of the function / with respect to the variable 
U . Such problems arise in population dynamics and Differential Geometry. The 
difficulty of studying this problem is that this nonlinear elliptic system does not 
fit the usual growth condition in M.Giaquinta's book [G] and the natural working 
space H 1 n L°°(0) for the corresponding Euler-Lagrange functional does not fit 
the usual minimization or variational argument. We use the direct method on 
a convex subset of H 1 n L°°(f2) to overcome these difficulties. Under a suitable 
assumption on the function /, we prove that there is at least one solution to 
this problem. We also give application of our result to the Dirichlet problem of 
harmonic maps into the standard sphere 



1. Introduction 

In the interesting paper, Bensoussan-Boccardo-Murat [BBM] studied a class of 
quasi-linear elliptic equations which include the following equation as special case: 

— Au + u\Du\ 2 = h(x), in Q u\gn = 0. (*) 

where Q a smooth bounded domain in R n with h G L°°()fl. Observe that equation 
(*) has a variational structure for left side of (*). In fact, let p(u) = —u 2 /2. Then 
the equation can be written as 

-e- p{u) div(e p{u) v«) = h(x), in Q, u\ dQ = 0. 

Let P(u) = J "e p(s) (is. Clearly P = P{u) has its derivative P'(u) = e p(s) > 0, and 
has its inverse function u = G(v). Define v = P{u). Then the equation can be 
reduced into 

-At; = e- G{v)2/2 h := h{x, v), in fi, v\ dn = 1. 
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In this case, the nonlinear term is in right side and it is very tamed. The variational 
structure is 



J{v) = \ [ \Dv\ 2 - [ H(x,v) 
A Jn Jn 



where H (x, v) = j" h(x, s)ds. Using the mountain pass method or variational tech- 
niques one can easily prove the existence of a classical positive solution of the prob- 
lem (*). This kind of idea is usually called the Cole-Hopf transformation in the 
study of parabolic equations. One may also obtain a Liouville type theorem for 
bounded smooth solutions of the problem: 

-e- p{u) div{e p(u) v u) = 0, in R n . 

If we consider the evolution equation of the following form: 

u t — Au + u\Du\ 2 = h(x), in Q, u\ t= o — ip, u\qq = 0. (**) 
We can also use the trick above and reduced this equation into 

vt - Av = e' G ^ 2 ' 2 h, inn, v\ t=0 = P(ip), v\ dn = 1. 

Hence one can also get a global existence of solutions of the problem (**). However, 
this trick is not so good for quasi-linear elliptic systems or parabolic systems. We 
can also consider (*) in another way. Let m(u) = —u 2 . Then the elliptic equation 
in (*) can be reduced into 

- e - m ^div(e m{u) v«) + \m'{u)\ v u\ 2 = h(x) 
and the left side is the variational derivative of the functional 



j( u ) = I e m{u) \Du\ 2 dx. 
Jn 



This tells us that we can use this functional to study the corresponding Dirichlet 
boundary value problem: 



- e - m ^div{e m{u) Xju) + -m'(u)\ XJ u\ 2 = 0, in Q u\ 
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Although the functional j(-) is not so nice, one can show that this problem is always 
solvable. 

In this paper, we mainly discuss the Dirichlet problem of the quasi-linear elliptic 

system 

- e -fMdiv{e f M \7U)+ l -f{U)\ V = 0, in SI, U\ dn = 4>. (1.1) 

Here again Vl a smooth bounded domain in R n , f : R N — > R is a smooth function, 
U : Q — > R N is the unknown vector-valued function, : dSl — > R N is a given 



vector-valued C 2 > 5 function, /' is the gradient of the function / with respect to the 

variable U . Such problems arise in population dynamics and Differential Geometry 

(see below). The difficulty of this problem is that this nonlinear elliptic system 

does not fit the usual growth condition in M.Giaquinta's book [G] (see also [C] 

and [S]) and the natural working space H 1 fl L°°(fi) for the corresponding Euler- 

Lagrange functional does not fit the variational argument. We point out that the 

space H 1 fl L°°(Cl) or W lA (q > 2) is the right for the corresponding variational 

functional to be differentiable. Once we have this, we can derive the Euler-Lagrange 

system for the functional. In our earlier research [M], we noticed that there is 

one nice direction for us to differentiate the variational functional. Then we use 

this property and the Nash-De Giorgi- Moser-Morrey iteration to obtain such a L°° 

bound. This method is tricky, we used it in the study of certain elliptic systems (see 

[L]). However, it is not so powerful. In fact, we can not use it to treat the system 

(1.1). Our idea here is to find a nice space obtaining a L°° bound. To get such a L°° 

bound, we use the direct method on a convex subset of H 1 nL°°(f2). In the following 

we will study the existence of a weak solution of (1.1) in the class H 1 fl L°°(Q). 

Theorem 1.1. Suppose that f(U) satisfies that f'(U) = —Ug(U), where g(U) is 
a positive continuous function on R N . Then there is at least one weak bounded 
solution U of (1.1). Furthermore, there is some q > 2 such that U G W 1,q (Q). 

Typical examples satisfying our assumption of Theorem 1 are (1) f(U) = — a\U\ 2 ,(2) 
f{U) = -21ogi±p, and (3) f(U) = -/31og±±^f where a, (3 > 0. Case (1) cor- 
responds to Gaussian measure dm = e~^ Y ^dY. In this case, f'(Y) = —2aY. Case 
(2) is for the standard metric ds 2 = ^ 1+ ^y^ i dY 2 on S 1 ^' 1 written in the stereo- 
graphic projection coordinates (Y) from the north pole. For this case, we have that 
f'(Y) = —AY /{I + \Y\ 2 ). So the critical map of E(-) is in fact a harmonic map 
from Q to the sphere. Then the result was obtained by S.Hildebrandt, H.Kaul, and 
K.O.Widman [HKW], Schoen and Uhlenbeck in [SU] and Giaquinta and Giusti in 
[GG]. When our / satisfies some convexity condition, one may consult the papers 
of S. Hildebrandt [H], D.G.Defigueiredo [D], and Marcellini-Sbordone, [MS]. As a 
by-product, we can obtain the following extension of a result due to V.Benci and 
J.M.Coron in dimension two [BC]. 

Corollary 1.2. . Given a bounded domain Q in R n with regular (Lipschitz) bound- 
ary. Assume we have a C 2,s map : dVL — > S N . Auume n < N. Then these exist 
at least two weakly harmonic maps from fl — > S N with Dirichlet boundary value 0. 

Proof. Since n < N and : dil — > S N is smooth, we have find a point P in S N such 
that {P, —P} are not in the range (j)(dQ). Taking P as the north pole and using the 
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stereographic coordinates at P, we can reduce the harmonic map problem into our 
problem (1.1) with f(U) = —2 log 1+ ^ . By Theorem 1.1, we have a weak bounded 
solution U, which corresponds a weakly harmonic map from Q — > S N with Dirichlet 
boundary value 0. In the same way, using — P as south pole, we can obtain another 
weakly harmonic map from Q — > S N with Dirichlet boundary value 0. □ 

Our method can be used to handle the following Dirichlet problem of the quasi- 
linear elliptic system 

-e- f ^d l {e f ^A{x)p 3 U a ) + l -[f {U)\ b A^d^U* = 0, in Q, U\ dn = 0. 

(l.iO 

Here we use the Einstein sum convention and we assume that (Af-) is a uniformly 
positive matrix function in the sense that its each component is non-negative and 
there are two positive constants A and A such that 

Mi? < A(x)$£f> < A^l 2 , (La) 
for any £ G R n x R N . The energy integral for this problem is 

E(U) = [ e f(u) A^diU a d 3 U b dx 
Jn 

Our result for this problem is 

Theorem 1.3. Suppose that f(U) satisfies that f'(U) = —Ug(U), where g(U) is 
a positive continuous function on R N . Assume that the matrix function (A^(x)) 
satisfies (l.a). Then there is at least one weak bounded solution U of (1-1 ')■ Fur- 
thermore, there is some q > 2 such that U G W l,q (Q). 

We also remark that we can extend our result to the Dirichlet boundary problem 
on half space i?" 

- e -fWdiv(eW v u) + !/'([/)! y u\2 = o, in ^ = . (L1 //) 

where is a smooth bounded function on i?". 
We have the following result: 

Theorem 1.4. Suppose that f(U) satisfies that f'(U) = —Ug(U), where g(U) is a 
positive continuous function on R N . Assume that is a smooth bounded function 
on i?" with f Rn \D<f)\ 2 dx < +oo. Then there is at least one weak bounded solution 

U of (1.1"). Furthermore, there is some q > 2 such that U G W^^R^). 

Clearly one can extend our results to the case when / = f(U) is replaced by a 
more general function / = f(x,U). One may also discuss the Dirichlet problem 
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of a quasi-linear elliptic system with a p-Laplacian type operator on the smooth 
bounded domain f2 in R n or in the half space it!": 

-L P U = in n, U\ dn = <p. (1.2) 

Here, p > 1, and 

L p u := e^^div^^DU^DU) + \f\U)\ y £/| 2 
with / e C 2 (R N ). Since the formulations are complicated, we omit it in this paper. 

2. Proofs of Theorem 1.1 and 1.3 

We first prove Theorem 1.1: 
We define the following integral 

E(U) = [ e^\DU\ 2 dx 
Jn 

Then one can formally compute that 

<DE(U), V >= [ e f{u \DU,DV) + - [ e f{u) f \U)\DU\ 2 dx 
Jn 2 J n 

Hence the Euler-Lagrange equation for E is (1.1). 

Since has a smooth extension over Q, we can find a constant vector C such that 

—C < (f>(x) < C for each x e Jl. Introduce the convex subset 

^ = {[/ e if 1 (Q); -C < U(x) < C, in n, U(x) = <f)(x), in dVl} 

Clearly <p E A and A is a weakly closed convex subset of H 1 . The functional E 
is weakly lower semi-continuous on A and C 1 on i/ 1 fl L°°(Q). Hence, we have a 
minimizer U E A such that E(U) = mi ve ^E(v). We will prove that this U is a 
weak solution of (1.1). Then we prove Theorem 1.1. 

Take any vector- valued function £ E Cq(Q) and a small positive constant e. Define 

V e (x) = min{C,max{—C, U(x) +e£(x)}}. 

Here the maximum and minimum are taken for each component. Then we have 

V £ E A. 

Define 

fi e := {x E Q; U(x) + e£(x) >C> U(x)} 

and 

ft e := {x E ft; C/(x) + e£(x) < -C < C/(x)} 
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Observe that for e > small we have U > on e and U < on Vt e . We also have 
the measure \Vt e \ — > and the measure f2 e — > as e — > 0+. 
Let 

f (x) = maa;{0, -C + C/(x) + ef (a)}. 

and 

£ e (x) = -min{0, C + U(x) + ef (a;)}. 

Clearly, £ e (x) > and f e (x) > 0, and both are in Hi n Then we can rewrite 

K as 

Since, for t e [0, 1), 

E(U + t(V e -U)) > E(U), 

we have 

< DE(U),V e -U>>0. 

Note that 

< DE(U), V e -U >=< DE(U), ef - f + £ e > . 
Hence we have that 

< DE(U),£ >> \[< DE(U),C >-< DE(U),Ce >] 

We will show that < DE(U),£ e >> o(e) and < DE{U),£ t >> o(e). In fact, we 
have, 

<DE{U),C> = [ e f{u) (DU,D£)dx 
Jn 

+ \ f(e^f'(U)-C\DU\ 2 dx 

= f e ^\DU\ 2 + - [ e^f(U)-(U-C)\DU\ 2 

+ ef e H u \DU,DZ) + le f f (U) ■ i\DU\ 2 

= / eW)\DU\ 2 g{U){l-\u.{U-C)) 

+ ef e^ u \DU,D()+ l -e f f (U) ■ £\DU\ 2 

> ef e W\DU,D0 + \e f e^f'(U) ■ t\DU\ 2 = o(e). 



In the last inequality we used the fact that for e > small we have U • (C — U) > 
on Vl e (since U > and C — U > on and in the last equality we used — > 0. 
Similarly we can prove that < DE(U),£ e >> o(e). Therefore, we have 

< DE(U),£ >> o(l) ->0+ 

Since our £ is arbitrary, we have 

<DE(U),£ >=o(l) ^0+ 

and 

>=0. 

Before finishing this section, we point out that our U is a spherical Q— minima in 
the sense M.Giaquinta and Giusti [GG]. Hence we have U G W 1,q (Q) for some q > 2. 
This proves our Theorem 1.1. 

The proof of Theorem 1.2 is a easy adaptation of the argument above. So we omit 
the detail. 



3. Proof of Theorem 1.4 
In this section we will give the proof of Theorem 1.3. 

We choose a bounded domain exhaustion of i?" in the way that i?" = £/^>o(-R™ ) D 
.Br(O). Here B R (0) is the ball of radius R with center at 0. We write = 

(^)nB fi (o). 

For any bounded domain Q in i?", we let 

£ n ([/) = I e m \DU\ 2 dx 
Jn 

Fix a large R > 0. We define the following integral on 

Er(U) = [ e fiu) \DU\ 2 dx 

By Theorem 1.1, we have a bounded weak solution ur with uniform bounds: 

I m rIl°°(_b+) ^ I m -r|l°°(r^) 

and 

I jDM kIl2(B+) < C 2 \D(f)\ L 2 (R n )} 

where C 2 is a unform constant. 



With these two bounds we can use extracting a diagonal subsequence method to 
get a weakly convergence sequence (U k ) and a weak solution W G -ff/ oc (i?")nL oc (i?™) 
of (1.1") with the following property 

E^W)<lmi k ^E^U k ). 

This proves Theorem 1.3. 
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